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Abstract. We present observations of the dwarf irregular galaxy NGC 4449 at 850 /.im and 450 jim obtained with SCUBA at 
the JCMT. The distribution of the cold dust agrees well with that of the CO and Ho- emission. To explain the integrated mm- 
through far-infrared continuum spectrum three dust components are required, with temperatures of 16 K, 39 K and 168 K, 
respectively. The dust mass is dominated by the cold dust component; we derive a total dust mass of ~ 3.8 x 10 6 M Q , and 
with the local gas-to-dust ratio of ~ 130 a total gas mass of M H i + h 2 ~ 4.9 x 10 8 M Q . Comparison with the HI mass leads to 
a total molecular gas mass of ~ 3.4 x 10 s M Q . We derive a conversion factor of the CO line intensity to molecular hydrogen 
column density X C o = Mfc/A:o which is at least 11 times larger than the Galactic value. These values are in accord with the 
lower metallicity of NGC 4449. 
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1. Introduction 

Carbon monoxide (CO) is commonly used as a tracer for 
molecular gas because the molecular hydrogen (H2), the dom- 
inant species, lacks strong emission lines from which the col- 
umn density could be determined. Because the rotational tran- 
sitions of CO are relatively easy to excite, it is possible to 
1 use the CO line intensity integrated over all velocities (Iqo = 
• J Tmbdv) or the CO luminosity, Leo, in one of these lines to es- 
timate the column density, NQI2), and mass, M(H2), of molec- 
ular hydrogen, provided one knows the correct conversion. For 
the Milky Way the conversion factor, Xco = Nh 2 /7co(1-0), has 
been found to be between Xco = 1.6 X 10 20 cm~ 2 (Kkms -1 ) -1 
(Hunter et al . [19971 and X C o = 2.3 x 10 20 enr 2 (K km s" 1 )" 1 
(Strong et al. 1988). Throughout the paper we will use the lat- 
ter because it is widely used for comparison in the literature. 
Early on, the question arose whether one can apply this local 
value to external galaxies, as the relation may strongly depend 
on the physical conditions such as the radiation field and the 
metallicity (Dettmar & Heithausen 1 19891 Israel 1 19971 Barone 
et al. 2000). Dwarf irregular galaxies provide excellent labo- 
ratories to investigate Xco in environments with lower metal 
abundances and varying excitation conditions. 

An alternative tracer for the molecular gas is the cold dust 
that one can measure at mm/submm wavelengths (Guelin et 
al. fT993l Dumke et al. fT997i . The dust mass can be deter- 
mined from the emission spectrum with some assumptions on 
the dust emissivity and composition. Adopting a gas-to-dust 
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ratio the dust mass can be used to estimate the molecular gas 
mass. This is however not a trivial business because the gas- 
to-dust ratio is derived locally in the Milky Way and may vary 
for different conditions. The CO abundance strongly depends 
on the amount of dust, because the lower the amount of dust, 
the lower the shielding of molecules, hence they are easier to 
dissociate, which in turn will render CO as an unreliable tracer 
of the molecular gas. 

To study the relation of cold dust and molecular gas 
in low metallicity environments we have therefore decided 
to observe NGC 4449 with the Submillimetre Common-User 
Bolometer Array (SCUBA) at the James Clerk Maxwell 
Telescope (JCMT) at 850 and 450 /mi. Maps of the cold dust 
in the millimetre and submillimetre regime for dwarf galaxies 
are rather rare. Together with the investigation of NCG 1569 
(Lisenfeld et al. I2002> the study of NGC 4449 presented here 
is the second of its kind. 

The galaxy NGC 4449 has seen numerous studies over es- 
sentially the whole electro-magnetic spectrum (see e.g. Hunter 
et al. rH%9l Hill et al. rPMl Klein et al. fT996l . It is a nearby 
Magellanic irregular dwarf galaxy of high surface brightness 
at a distance of about 3.7 Mpc (Bajaja et al. 119941 . The recent 
investigation of Karachentsev et al. (2003) yielded a distance 
of 4.12+0.50 Mpc, which agrees with the older value within 
the errors. Located in the constellation Canes Venatici it is a 
member of the CVnl group (Sandage I1975> . which contains 
mainly spirals and irregular galaxies. The observed morphol- 
ogy and kinematics of the HI gas suggests a tidal interaction 
with the nearby dwarf galaxy DDO 125 (Theis & Kohle l200Tl 
and may explain the high star-formation rate and the huge HI 
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streamers surrounding NGC 4449. The metallicity was found 
to be 12+log(Q/H) = 8.3 (Lequeux et al.. lT979l 

Hunter et al. (1986) stated that the spectral distribution 
of the far-infrared emission from NGC 4449 as observed with 
IRAS is not different from normal spirals. The ratio of the 
IRAS flux densities at 100 and 60/mi S ioo/jm/S60//m ~ 2 is at 
the lower end of the range covered by spirals and indicates a 
high star formation rate. However, the dust mass calculated 
from this measurement only relies on the warm dust compo- 
nents, thus representing a lower limit. Our new observations 
allow to determine the amount of cold dust in NGC 4449 and 
to study the gas-to-dust ratio in more detail. In Sect. 2 we de- 
scribe the observations and data reduction procedure. Sect. 3 
presents the results, while in Sect. 4 a comparison with the gas 
in NGC 4449 is made. A summary is given in Sect. 5. 

2. Observations and data reduction 

NGC 4449 was imaged simultaneously at 850 and 450 fim us- 
ing the SCUBA camera at the JCMT in January 18, 1999. 
The observations were made in the jiggle mode with a chop- 
per throw of 120". Two fields of view of about 23 with suffi- 
cient overlap were observed. The extent of the galaxy at optical 
wavelengths is about 2' by 3'. Since we observed about 9 hours 
with the source rotating on the sky, there may have been some 
emission from the galaxy itself in the reference beam occasion- 
ally. However, by combining all coverages this error should be 
small compared to the calibration uncertainties. Thus only a 
diffuse extended dust component possibly associated with the 
large HI halo would have been missed by our observations. 

The total net observing time was 8960 s for field A and 
6528 s for field B. The two fields were placed around a central 
position atai9 5 o=12 /, 25'"46 s , (5i 950 =44 o 22'45" such as to cover 
the main (CO emitting) body of the galaxy. The whole observ- 
ing run (9 hours in total) had good and stable weather with an 
optical depth of t^q ss 0.23 and a seeing of ~ 074. Flux calibra- 
tion was made using Mars, which was observed in the middle 
and at the end of the observing run. 

The data have been reduced in the standard manner using 
the STARLINK reduction package SURF (Jenness & Lightfoot 
1999 1. Flatfielding was performed to correct for differences in 
the sensitivity of the bolometers. The atmospheric opacity was 
determined from skydip observations. Noisy bolometers and 
spikes were identified and blanked out. Typical pointing cor- 
rections were 1" to 2", hence small compared to the beam size. 
From the planet profile, the beam size (HPBW) was estimated 
to 1573 and 875 at 850 and 450 fim, respectively. Fig.^shows 
the resulting maps of NGC 4449 at 850 and 450 fim in the top 
panels. The noise levels of the maps are <x = 6 mJy/beam at 
850 fim and cr — 18 mJy/beam at 450 fim. 

Observations at 1200 fim were previously carried out with 
the IRAM 30 m telescope on Pico Veleta (Spain) in December 
1998, using the 19-channel bolometer array of the Max-Planck- 
Institut fiir Radioastronomie (Kohle 1999 1. We obtained On- 
The-Fly (OTF) maps with sizes between 300"x200" and 
400"x300" in azimuth and elevation, with the telescope mov- 
ing in azimuth at a speed of 6" /s. Successive horizontal scans 
were spaced by 4" in elevation. The wobbler throw was 45" 



and wobbled at 0.5 Hz in azimuth. During the run the weather 
conditions were good and the atmosphere stable, with sky 
opacities between 0.1 and 0.3 at 230 GHz. Two channels had 
to be discarded. The total integration time was about 8 hours. 
Calibration was made using Mars, which was observed twice 
during the run. The derived beam-size is 1078 and the resulting 
map is shown on the left of the centre row in Fig.^ smoothed 
to match the 1573 of the 850 fim map and to improve the signal- 
to-noise. The noise level of this map is cr = 2 mJy/beam. 

We have to mention that restoring maps from chopped on- 
the-fly measurements may cause photometric imperfections, as 
the information on the large-scale structure (larger than the 
wobbler throw) is lost. The same holds true of course for the 
IRAM data: we may have missed a diffuse, extended dust com- 
ponent and hence a noticeable fraction of the integrated flux. 
Such a component was suggested by Thronson et al. ( 1987 > 
which detected emission at 150 fim with the Kuiper Airborne 
Observatory (KAO) on scales of up to 4 arcmin at a low signal- 
to-noise ratio. However, this will not change our general con- 
clusions. In particular, if we would be able to correct for a miss- 
ing large-scale component, we would flatten the slope of the 
submm spectrum, and increase the dust mass (see Sect. 3). 

3. Results 

Fig. 1 exhibits the dust continuum emission of NGC 4449 re- 
sulting from our observations, along with images in the HI, Ha 
and CO line. The galaxy is clearly detected at all mm/submm 
wavelengths. The morphology of the galaxy is rather similar for 
all components, except for the neutral atomic hydrogen. The 
strongest dust emission comes from a bar-like structure promi- 
nent in the optical regime, although the brightest regions do not 
coincide. The optically bright northern region is conspicuous in 
all maps and manifests a large star-forming complex which is 
most pronounced in the Ha map. A more detailed comparison 
is hardly possible, due to the low signal-to-noise and therefore 
the uncertainties in the exact structures. 

By integrating over the whole SCUBA maps we have de- 
termined total flux densities to be Ssso/jm = 1-8 ± 0.2 Jy and 
= 16.5 ±4.1 Jy. No color correction has been made. 
Because the CO(3-2)-line emission falls into the bolome- 
ter bandpass at 850 /urn, we corrected the total continuum 
flux for the line contribution. We estimated the contamina- 
tion using the CO(1-0) map and assuming a line ratio of 
^co(3-2)/^co(i-0) ~ 0.7 (Fritz 2002, private communication). 
We thus obtained a contribution of the CO line to the to- 
tal emission of ~ 12%, resulting in a pure continuum flux of 
1.5 ± 0.3 Jy. The contamination at 450 fim is negligible com- 
pared to the strong thermal dust emission, hence no correction 
was made. The integrated flux at 1 .2 mm was taken from Kohle 
( 199 < 5Jl and is also corrected for contributions of the CO(2-l) 
line. The quoted errors are formal errors based on the map noise 
and on the uncertainty of integrating the flux, including calibra- 
tion errors. 

In the literature there are 5 published FIR flux densities 
for NGC 4449. These are the four IRAS data points between 
12 fim an 100 fim (Hunter et al. 1986 1, and a 150 fim value 
from KAO observations (Thronson et al. I1987l >. Together with 
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Fig. 1. Intensity maps of NGC 4449 at 850 (upper left) and 450 fim (upper right), presented in this work, with contour levels 
of 2,3,4 and 5 cr, where cr is the noise level of the map (cr(450/mi) = 18mJy/beam and cr(850jum) = 6mJy/beam); at 1.2 mm 
(centre left, from Kohle 1999), with contour levels of 2, 2.5 and 3 cr (cr(1200jum) = 2mJy/beam); an Ha image (centre right, 
from Bomans et al. 1997 1, the integrated HI (bottom left, from Hunter 1999 1, with contour levels of 0.06, 0.08 and 0.1 Jy/beam 
and the integrated CO(1-0) (bottom right, from Kohle et al. 1998 1999), with contour levels from 0.5 to 2.5 K km s _1 in steps of 
0.5 K km s _1 . All maps show the same area as that mapped with SCUBA. 
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Table 1. IR to submm integrated flux densities of NGC 4449 
from 12 jum to 1200 fim. The values at 850 and 1200 fim have 
been corrected for contamination by CO lines. 
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A 


flux density 






(Jy) 




12 


2.1 ±0.2 


Hunter et al. (1986 1 


25 


4.7 ± 0.5 


Hunter et al. fT986l 


60 


36 ±3 


Hunter et al. <1986> 


100 


73 ± 8 


Hunter et al. (1985) 


150 


100 ± 20 


Thronson et al. RW!\ 


450 


16.5 ±4.1 


this work 


850 


1.5 ±0.3 * 


this work 


1200 


0.26 ± 0.04 * 


Kohle (19991 



* corrected for CO line contamination 
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Fig. 2. The mid-infrared to millimeter spectrum of NGC 4449, 
fitted with three-temperature dust components. Crosses mark 
flux densities from Table 1. The three dust components have 
temperatures of 168 K (dashed line), 39 K (dotted line) and 
16 K (dashed-dotted line). 



the data points presented in this paper this leads to a fair sam- 
pling of the spectral energy distribution of NGC 4449 from mm 
through FIR wavelengths. In Table 1 we have compiled all val- 
ues, Fig. |2 showing the resulting spectrum. 

If we assume the dust grains to be in equilibrium with the 
radiation field we can describe their emission by a modified 
black-body law. This will definitively be the case for large 
grains, which contribute most to the sub-millimetric contin- 
uum. They also account for more than 90% of the dust mass 
(see e.g. Desert et al. [T990l. The emission from 20 to 80 fim 
can originate from very small grains at a significant amount. 
In this case, with our simple approach we may overestimate 
the mass of that component by a few percent. However, as we 
will see below, the mass of this hot component is of almost no 
importance. 

In case of a low optical depth the intensity is given by 
F A = B A (T) ■ k(A), where B A (T) is the Planck function, T the 
dust temperature and k(A) is the absorption coefficient. The 
wavelength dependence of the latter, which is usually assumed 



to follow a power-law k(A) oc A~P (see e.g. Alton et al. 2000 1, 
is not precisely known. Models of the interstellar dust yield a 
spectral index B — 1 ... 2 for mixtures of small and big grains 
(e.g. Draine & Lee 19841. Gordon ( 1988) has noted a tempera- 
ture dependence of B, for the range 1 ... 2.5, and Dunne & Eales 
( 2001 ) have derived B =s 2 on a statistical basis for a sample of 
bright IRAS galaxies. 

In order to further analyze the dust components of 
NGC 4449 we have modeled the spectrum with three dust com- 
ponents at different temperatures. This yields a reasonable fit 
to the emission spectrum over the entire observed wavelength 
range. As we have only eight observed spectral data points, we 
have to limit the model to three components. In reality, the dust 
is probably at a range of temperatures, reflecting the range of 
colours and intensities of the heating radiation fields. We de- 
cided to use B — 1 for the hot and B — 2 for the warm compo- 
nent, leaving 7 free parameters for the fit. It was not possible 
to fit the steep slope of the spectrum satisfactorily with B — 2 
for the cold component, thus we decided to leave it free up to 
2.5, resulting in a value of B — 2.3. Fig. 2 shows the result with 
hot (T = 168 + 10 K,f3 = 1), warm (T = 39 ± 3 K,B = 2) and 
cold (T = 16 + 2 K, B — 2.3) dust components. 

The derived parameters are in accord with other investiga- 
tions (e.g. Dumke et al. I1997L and agree well with measure- 
ments of the diffuse ISM in the Galaxy (T * 17 K,B = 2) by 
COBE/FIRAS (Reach et al. ll995l Sodroski et al. [TW7t . In or- 
der to properly fit the relatively strong (sub)millimeter emis- 
sion, a cold component is required. Moreover, because of the 
low dust emissivity at this temperature, a very large fraction of 
the total dust has to be cold in order to explain the emission 
at 850 fim and 1200 /mi. An attempt to fit a two-temperature 
model yields significantly worse results. 

From the three-component fit to the spectrum of NGC 4449 
we can estimate the dust mass. This can be achieved with a 
simple equation adapted from Hildebrand ( 1983 1, where 

D 2 Ss50fi m 

Md = b ■ (1) 

«850/jmU )*850/fln 

Here, B^o^miT) is the Planck function at 850 fim, D is 
the distance and S l 850 J um the observed flux density. To derive 
a consistent value of the dust absorption coefficient k^o^ we 
used the above mentioned power-law k(A) oc A~P and the canon- 
ical value of feo/jm = 10 cm 2 g _1 from Hildebrand ( 1983), ex- 
trapolating to the relevant wavelengths and using the asso- 
ciated value of B. We thus obtained the following values of 
the dust absorption coefficient for the three different compo- 
nents: = 0.6 cm 2 g _I for the cold (16 K) component 
at 850 fim, fcioo//m = 63 cm 2 g -1 for the warm (39 K) and 
&ioo/im = 25 cm 2 g -1 for the hot component at 100 fim. These 
values are consistent given our assumptions, and especially for 
the spectral indices used here. However, in the literature one 
can find a wide range of values for the emissivity and the de- 
pendences are only poorly known. As the most extreme cases 
we quote two values, namely k^o^ = 0.26 cm 2 g -1 (Alton et 
aL EDOa and k S5Qflm = 2.4 cm 2 g" 1 (Alton et al. 120001 . 

The three components shown in Fig. |2J result in a very 
small mass of hot dust, M^h -18 M Q , a larger mass of warm 
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dust, Md jW = 0.3 x 10 5 M and a large mass of the cold dust, 
M d>c = 3.8 x 10 6 M . The total dust mass, M d = 3.8 x 10 6 M G , 
is thus completely dominated by the cold dust component. The 
uncertainty of this mass is naturally large (within a factor of 
3), and is mainly governed by the peak and the slope of the 
spectrum between 100 - 1200 yum, and by the assumed dust ab- 
sorption coefficient, which is not well known. The error in cal- 
culating the dust mass caused by using modified Plankians is 
negligible in view of the large uncertainties in the general grain 
properties. 

Thronson et al. (Q987 ) calculated a somewhat lower mass 
of M d a; 1.5 x 10 6 M , using their 150 //m value and assuming 
a spectral index of /3 — 1 ; the deviation to our value is most 
likely caused by the then unknown slope of the spectrum in the 
millimetre regime. 

Lisenfeld et al. (2002i have undertaken a similar study of 
NGC 1569, a galaxy which is presently in the aftermath of 
a massive burst of star formation. The dust spectrum of that 
galaxy shows a completely different slope in the range between 
150 - 1200 /mi. The IRAS-flux densities of NGC 1569 indicate 
a different location of the maximum between 100 and 60 //m, 
compared to NGC 4449, indicating the presence of much more 
hot dust in the former. A very cold component of around 7 K 
was needed to properly fit the spectrum, which resulted in a 
tremendous mass of this very cold dust. The authors stated 
that such a large amount of cold dust is very unlikely in a 
dwarf galaxy dominated by intense radiation fields (2 orders 
of magnitude higher than in the solar neighbourhood). Hence, 
Lisenfeld et al. favour the model of Desert et al. ( 1990 1, where 
the very small grains have optical properties consistent with 
beta=l. Since NGC 1569 has a shallow slope, the resulting fit 
requires a large contribution of grains emitting like beta=l, 
thus a large contribution of VSG with the Desert's et al. prop- 
erties. The dust mass derived from this model is much less 
than that from a fit with modified Plankians and results in a 
overall gas-to-dust ratio of M g /M d w 1500 ••• 2900. We have 
to point out, however, that the physical reason for the steeper 
slope of the submm dust continuum spectrum of NGC 4449 
08 = 2.3 vs. f3 = 1.0 in case of NGC 1569) is not at all ob- 
vious. 

4. Comparison with the gas 

The derived dust distribution can be compared with other trac- 
ers of the ISM, viz. CO and HI. We also have a CO(1-0) map 
of NGC 4449 at our disposal, obtained with the IRAM 30-m 
telescope on Pico Veleta (Kohle ll998lfT999b . Fig.Q]shows the 
integrated map, the angular resolution is 22". The correlation 
of CO and dust is clearly seen. However, the CO coverage of 
the galaxy is not uniform, especially in the outer parts so that 
not all of the structures seen there are reliable. 

We have also compared the Ha map published by Bomans 
et al. J1997I and shown in Fig. ^ to the maps of the dust and 
the molecular gas. It is evident that the strongest dust emission 
follows the distribution of intense star formation (as expected). 
However, this is not a simple one-to-one relation. The same is 
true for the CO emission and may indicate dissociation pro- 
cesses as observed in massive star forming regions. There is 



also significant dust emission away from the star-forming re- 
gions, but this is accompanied by CO radiation and may show 
a more diffuse component (smaller clouds) of the molecular 
gas. 

A comparison with the HI (Hunter et al.Q999l shows only 
little correlation. The inner HI of NGC 4449 is organized in a 
ring-like feature around the optical body. In the northeast one 
can see part of this ring, whereas in the south there is nearly no 
HI within our mapsize. In the center of NGC 4449 only little HI 
is found, which is expected. The gas should be ionized there, 
in molecular form or frozen out onto the dust grains. 

With the dust mass of M d = 3.8 x 10 6 M derived above we 
can estimate the total gas mass. Using the normal gas-to-dust 
ratio of M g /M d « 130 (Devereux & Young Ft^O! this yields a 
total gas mass of Mhi+h, = 4.9 x 10 8 M for the inner (optical) 
part of the galaxy. Computing the HI mass within our map- 
size yields Mui = 1.5 X 10 8 M . Subtraction from the total gas 
mass estimated above results in a mass of molecular hydro- 
gen of Mh, = 3.4 x 10 8 M . From the integrated CO spectrum 
we derive a CO luminosity of Leo = 7.9 x 10 6 Kkm s _1 pc 2 , 
which leads to a conversion factor from CO intensity to H2 col- 
umn density Xqo= 25 x 10 20 crrT 2 (Kkm s -1 ) -1 , i.e. about 11 
times the standard Galactic value quoted in Sect. 1. Such an 
enhanced value is explicable in view of the lower metallicity, 
of NGC 4449, viz. 12+log(0/H) = 8.3 (Lequeux et al.. 179791 . 

Kohle (1999 1 derived an even higher Xqo, assuming 
virialized cloud complexes in the CO map. He obtained 
X co = (37 ± 11) X 10 20 enr 2 (Kkm s -1 ) -1 , which is about 16 
times the Galactic value and may represent an upper limit. 
This would lead to a much higher molecular gas mass and 
hence to an even higher gas-to-dust ratio. Both scenarios are 
in agreement with what has been found for other low-mass 
dwarf irregular galaxies (see e.g. Barone et al., 2000 and refer- 
ences therein). Thronson et al. ( 1987) derived Xco /^mw ~ 10 
from their KAO data. Based on a much poorer sampling of the 
dust spectrum using IRAS CPC data Israel ( 1997) calculated 
a value of Xco = 8 x 10 20 crrT 2 (K km s ) _1 , which is signifi- 
cantly lower than the other values discussed here. 

5. Conclusions 

We have presented new maps of the dwarf galaxy NGC 4449 
at 450 and 850 jum, taken with SCUBA at the JCMT The in- 
tegrated flux densities were compared to those obtained with 
IRAS and KAO. The derived flux densities can be readily ex- 
plained by the presence of a significant amount of cold dust. We 
were able to fit the observed SED well with modified Plankians 
of three dust components, requiring most of the dust to be at a 
temperature of 16 K, with a spectral index /3 = 2.3, which is a 
bit higher than the frequently used value of 2. 

Comparison with CO observations shows a good coinci- 
dence of the cold dust and the molecular gas. The Ha image 
shows that star-forming regions are related to the dust emission 
peaks. The main emission comes from the bar-like structure, 
but there is also weak dust emission away from the main body 
of NGC 4449, and there may be some more following the HI 
ring-structure. 
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Assuming a standard local gas-to-dust ratio, we derive a 
conversion factor Xco ~ 25 x 10 20 cm -2 (Kkm s -1 ) -1 , which 
is 11 times the Galactic value. Given the low metallicity of 
NGC 4449 (1/3 solar), this is probably still a lower limit, be- 
cause a standard gas-to-dust ratio may not be applicable. The 
standard Galactic conversion factor thus strongly underesti- 
mates the molecular gas mass in such galaxies, because our 
tracer - the CO molecule - is less abundant, due to the lower 
metallicity and/or dissociation. In view of the low metallicity 
of NGC 4449 this seems in perfect agreement with studies of 
other dwarf irregulars. 
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